A GAUGE THEORETIC APPROACH TO ELASTICITY WITH 

MICROROTATIONS 



C. G. BOHMER AND YU.N. OBUKHOV 



Abstract. We formulate elasticity theory with microrotations using the frame- 
work of gauge theories, which has been developed and successfully applied 
in various areas of gravitation and cosmology. Following this approach, we 
demonstrate the existence of particlc-likc solutions. Mathematically this is 
due to the fact that our equations of motion are of Sine-Gordon type and thus 
have soliton type solutions. Similar to Skyrmions and Kinks in classical field 
theory, we can show explicitly that these solutions have a topological origin. 



1. Introduction 

1.1. Elasticity and gauge theories. In classical elasticity one assumes that the 
elements of a material continuum are structureless points which are characterized 
by their positions only, whereas their internal properties (such as orientation, defor- 
mation etc) are disregarded. This picture is generalized to the concept of continua 
with microstructure when a material body or fluid is formed of a large number of 
deformable particles whose (micro)elastic properties contribute to the macroscopic 
dynamics of such a material medium. The corresponding mechanical model was 
developed very early by the Cosserat brothers [8] who replaced a structureless ma- 
terial point by a point that carries a triad of axes (an oriented frame, in the modern 
geometrical language) . This idea has led to a rich variety of models which are known 
by various different names like oriented or multipolar medium, asymmetric elastic- 
ity micropolar elasticity see e.g. [ll[T5l[ll[17l[2ll|43i55,^^^^tl9j. This 
field is also related to other areas of continuum mechanics, like the theory of fer- 
romagnetic materials, cracked media, liquid crystals 0, superfluids [65 l l46 l [58 l 166] 
and granular media (see, for example, application of the Cosserat model to the 
description of the grain size effects in metallic polycrystal [68j). 

In general case, the 3-frame attached to material elements is deformable, thus 
introducing nine new degrees of freedom in addition to the classical elasticity (3 
rotations, 1 volume expansion, and 5 shear deformations); this is usually called a mi- 
cromorphic continuum [19j . A special case is the microstretch continuum for which 
the 3-frame can rotate and change volume (4 degrees of freedom: 3 microrotations 
and 1 microvolume expansion). Finally, in the micropolar continuum, the struc- 
tured material point is rigid, thus possessing 3 microrotation degrees of freedom. 
The last is the subject of the Cosserat elasticity that considers media which can ex- 
perience displacements and microrotations. In turn, this theory has two important 
limiting cases, one with no microrotations which is the classical elasticity, and the 
second case where one assumes that the medium only experiences microrotations 
and no displacements. Models of the latter type have in fact a long history which 



2000 Mathematics Subject Classification, (primary), (secondary). 

1 



2 



C. G. BOHMER AND YU.N. OBUKHOV 



can be traced back to MacCuUagh in 1839, see [67l|64]- Recently OH] such models 
have been investigated in the fully nonlinear setting and plane wave type solutions 
were explicitly constructed. The existence of such solutions is a highly non-trivial 
fact. In a linearised setting, similar solutions were investigated in [48l l49l [50] . 

Recent years have seen a revival of elasticity, in particular in the theory of 
dislocations which has been analysed from a gauge theoretic point of view. The 
gauge approach was developed and successfully applied in high energy physics and 
in the theories of gravity [531 113 HH] • The gauge approach in elasticity also has a 
long and rich history, see for example [Ml[Illllll3[31[31[31IMl[3ilinillIl- Of 
particular interest is the fact that the dislocations behave very much like particles. 
From a more formal point of view, one of the open questions in this field of research 
is whether or not it is possible to find soliton type solutions which could justify the 
particle interpretation from a mathematical point of view. The main result of this 
paper is that wc are able to find soliton type solutions in a particular model of 
elasticity motivated by gauge theories of gravity. 

1.2. Geometry and elasticity. Our notation and conventions follow the book 
|26j where the modern differential geometry of manifolds is presented. Here we 
briefly summarize the notions and define objects needed for our current study. 

We consider an elastic medium which occupies the whole of K'^, and we can 
identify material points with points in space. We use the Latin alphabet to label the 
coordinate indices, i,j, fc, • • • = 1, 2, 3, say, we write for the spatial (holonomic) 
coordinates of A^, a 3-dimcnsional simply connected manifold which we identify 
with the material points of the continuum. This manifold occupies the entire 3- 
dimcnsional Euclidean space M."^. Greek indices refer to the (co)frame indices (also 
tangent space or anholonomic indices), we use a, /3, • • • = 1, 2, 3 to label the frame 
covectors (1 forms) ~ hfdx^ which we can refer to as the distortion in elasticity. 
We follow the Einstein summation convention whereby we sum over twice repeated 
indices. We denote the frame by Ca = h^^d/dx^, such that ea\^^ = h^^h^ = S^. 
The Ca form the basis vectors of the tangent manifold TAi, and J denotes the 
interior product. 

The Euclidean metric g is defined on as a scalar product g : TA4 x TAi — > M. 
As usual, the values of this scalar product on the basis vectors introduce the metric 
tensor gij = dj), the components of which for the Cartesian coordinates read 
Qij = Sij = diag(-|-l, +1, +1). The components of the metric gap — gifia^ S/s) in an 
arbitrary basis are obviously related to the frame components via gap = haSijh-'^. 
One can think of i?" as a set of three mutually orthogonal and normal covectors 
which form a basis at any point p € T*Ai. Each such basis covector has components 
(i?")i and thus we can also view the object hf as an orthogonal matrix. We can 
think of the metric of the deformed medium as the Cauchy-Green tensor in elasticity. 
Technically, the metric is used to raise and lower the indices, for example we obtain 
a covector- valued 1-form from a vector- valued coframe as -da = Sap'^^^ etc. The 
components of the coframe comprise the 3x3 matrix /if which is inverse to matrix 
composed of the frame components /i^ . From the definitions we directly obtain the 
following identity hf = g°'l^ Sijh-'p. 

One can draw a close analogy between gravity and elasticity, sec the instructive 
discussion [3S], for example. The force and the hyperforce (or the couple-force), 
as well as the stress (force per area) and the hyperstress (couple per area), caused 
by the force and hyperforce, respectively, are the basic concepts in the continuum 
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mechanics of structured media. The direct gravitational counterparts of the stress 
and hyperstress are the energy-momentum current and the hypermomentum (that 
inchides spin, dilation and shear currents). The spacetime itself is geometrically 
modelled as a continuum deformed under the action of the sources with nontriv- 
ial "stress" (energy-momentum) and "hyperstress" (spin current, e.g.). Using this 
correspondence, it is possible to put both gravity theory and elasticity theory into 
the same differential-geometric framework in which the metric, coframc, and con- 
nection are the fundamental field variables. We refer the reader for more formalism 
and technical detail to the comprehensive review |28] . 

When considering deformations of an elastic material, one places certain compat- 
ibility conditions on the induced stress. These so called Saint- Venant compatibility 
conditions are a form of integrability conditions so that the stress can be expressed 
and the symmetric derivative of the displacement. In a more geometrical language, 
by deforming the medium, wc do not want to induce any curvature into the mate- 
rial. This means that one assumes the Riemann curvature tensor of the deformed 
medium to vanish identically. It is well known that these two views arc equiva- 
lent, see for instance [TJ[60]. In the continuum theory of media with defects [31], 
the elastic properties of a crystal are described by means of the non-Riemannian 
geometries in which the notion of curvature is supplemented by a new geometrical 
quantity called torsion which was first introduced by Cartan in 1922. Therefore, we 
will analyse elastic materials using the language of the Riemann-Cartan geometries 
with torsion and we will stick to the well established requirement of having vanish- 
ing total curvature. The latter condition is sometimes called a second basic law in 
the differential geometry of media with defects |31| , and it has the physical meaning 
that the crystal structure is uniquely defined everywhere, see also pi l35[[5S] . In the 
following subsection we will collect the most important facts of such geometries. 

1.3. Geometries with vanishing Riemann curvature tensor. Let us describe 
each material point of our elastic material by a coframe i?" (a vector- valued 1-form) 
and a connection (a matrix- valued 1-form). Recall that the latter arises from 
the computation of the covariant derivative Vt, along a vector field v of the frame 
vectors, V„e/3 = r^"(w)eQ, see Sec. C.l of [26]. The connection is assumed to 
be metric-compatible which means the vanishing nonmetricity Dgap = 0, hence 
the connection 1-form is skew-symmetric, Fq,^ = — F^q (where F^^g = Ta^ gp~^). 
The coframe specifies the orientation of the orthonormal basis vectors at this point, 
while the connection determines how an arbitrary vector is parallelly transported 
near this point. The vanishing of the Riemann-Cartan curvature tensor (a matrix 
valued 2-form) takes the form 

(1.1) Rc,'^ ■.= dTj +tJ ST^" = Q, 

where d denotes the exterior derivative and A is the exterior (alternating) product. 
In geometries where the Riemann-Cartan curvature identically vanishes the con- 
nection is often referred to as a Weitzenbock connection. The condition Ra^ = 
is known in the literature by various names like teleparallel, fernparallcl of distant 
parallel geometry. Geometrically this condition means that the notion of paral- 
lelism is no longer a local statement, but can be defined globally, this means on 
the entire manifold. A straight line connecting two material points defines an abso- 
lute notion of parallelism. Physically the condition ()1.1|) has the meaning that the 
structure of an elastic medium is uniquely defined everywhere in M. , [35l IMl [31] [2] . 
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The only nonvanishing geometrical object in such geometries is Cartan's torsion (a 
vector valued 2-form) defined by 



(1.2) 



Let e ~ i)^ A i!)'^ A i!)^ denote the volume 3-form, then the dual forms of the 
products of the coframe are defined by 



(1.3) 
(1.4) 
(1.5) 



*(i?Q A dp A 7?^) = e^\eai3- 



The operator *, introduced by the formulas above, is called the Hodge operator. The 
dual forms satisfy the following useful identities 



(1.6) 
(1.7) 
(1.8) 



Aea 
d^ A e^. 
A Sctfii/ 



The object ScyPi is the totally antisymmetric Levi-Civita symbol with £123 = 1. 



1.4. Elastic invariants. Since is the only nontrivial geometric object, it is 
therefore the only non-trivial object characterising the deformations of the medium. 
It is natural to decompose it into its irreducible pieces which will serve as our 
building blocks of elastic invariants. The torsion T" = ^^"a"' Ad^ is represented 
as a rank 3 tensor which is skew-symmetric in one pair of indices, T^^j/ = —T"'^^. It 
has 9 independent components. Since we work in M? , it is also natural to consider 
the Hodge dual of torsion which can be regarded as a 3 x 3 matrix with no a 
priori symmetries. At any point we can view a 3 x 3 matrix as an element of a 9 
dimensional vector space V . Performing rotations which leave the metric invariant 
will change the components of this matrix, thus the group S0(3) acts on V and 
we are interested in identifying invariant subspaces. Those give rise to the so- 
called irreducible pieces of either this matrix, or its Hodge dual, the torsion tensor. 
Therefore, let us define T" := 

The three irreducible pieces of torsion are defined by 

(2>r" = i^"A(e^jr^), 

(3)r" ^ 1 *{dpAT^)e'', 
o 

As we already explained, the metric is used for the raising and lowering of the 
indices, i.e., e" = g°'^sp and -dp ~ gp-yi^^ ■ 

The irreducible pieces (|1.9|) - (|1.11|) can be conveniently understood in tensor lan- 
guage as the skew-symmetric part, the trace part, and the traceless symmetric 
part of the matrix, see the discussion below and in particular the explicit formulas 
(12:811- (I2T011 . 

It is important to note that there is a gauge freedom in choosing the variables 
Fq'^}. All coframes are related to each other by means of local rotations. 



(1.9) 

(1.10) 
(1.11) 
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Therefore, the same elastic medium can be described equivalently by another pair 
{•&'", T'^^} which is related to the former by means of the gauge transformation 

(1.12) = L"^i9'", 

(1.13) r;j" = L", r/ {L-^Yf, + l", d{L-^Yp, 

where L" ^ is an arbitrary 3x3 rotation matrix, this means L G S0(3). The 
covariant condition R^^ = is preserved by this gauge transformations, whereas 
torsion transforms covariantly, T'" — L'^a T"^- 

Making use of the local transformations (|1.12p - (|1.13p . we can eliminate either 
of the variables in the pair {d^jTa^}. The ability to do this is closely to our 
assumption of a globally flat manifold (jl.ip . In particular, in [3], the teleparallel 
connection is gauged away Ta^ = 0, and the coframe i?" is left as the only variable. 
Alternatively in [4] the coframe is gauged to the standard constant values hf = 6" 
and the only variable is the connection, in which case once can write 

(1.14) r^"=A",rf(A"i)^^. 

Here A G SO (3) parametrises a teleparallel connection. We will refer to this choice 
of variables as the connection gauge since the only dynamical variable is the con- 
nection. 

2. Basic equations of the model 

2.1. The potential energy. The natural potential energy is based on the sum of 
the three quadratic microrotational elastic invariants constructed from the nonva- 
nishing torsion 

3 

(2.1) y = ^ Q('>r"A*T„, 

1=1 

where c^, i = 1, 2, 3, are non-negative constants which are physically interpreted as 
elastic moduli. The non-negativity of the material constants Ci is a consequence of 
the condition of the positive semi-definiteness of the microrotational clastic energy, 
> 0, in complete agreement with [39] . 

In general, V — V{T"-) can depend arbitrarily on the torsion, but the choice of 
the quadratic function ()2.1|) yields the linear constitutive law, [39l|40l|4T]. When 
our medium occupies the whole of M'^, the three quadratic terms in p.l|) are not 
independent because of the teleparallel condition (|l.ip . Recall that there is an 
identity (see Eq. (5.9.18) of [55]) that is valid in n-dimensional space: 

(2.2) 

A Eap = R"'^ A - [(i>r" - (n - 2)(2>r" - i (3)j.a] /^^T^_d (2z?" A *r„) . 

The tilde denotes the Riemannian geometric objects, i.e. those constructed from 
the Christoffel symbols of the corresponding Riemannian metric of the manifold. 

Taking into account that the Riemann-Cartan curvature vanishes due to the 
teleparallel constraint (jl.ip and that by assumption the elastic medium is embedded 
in a flat Euclidean space with i?"^ =0, for n = 3 the identity p.2p allows us to 
express the square of the first (tensor) irreducible part in terms of the trace and 
axial trace squares: 

(2.3) (^>r" A *T„ = ^ ^ 1 (3)j.a ;^^To,~d (2i?" A *r„) . 
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In tensor language, using (|1.9p . (|1.10|) . and (|1.6p . (|1.7p . we have 

(2.4) (^ypa^^j,^ ^ ly.^^y^oMg^ 

(2.5) A = ^ (T^p.e^"")' e. 

It is straightforward to estabhsh relation between the irreducible parts of the 
torsion 2-form T" and its dual 1-form T". We have T" = T°'^'dp, denote the 
matrix T"^ by T, and we find explicitly the decomposition of the second rank 
tensor 

(2 7) T = '■"'"■'T + '•^^T + '•'^^ 

into the trace, antisymmetric and tracclcss symmetric parts: 

(2.8) (3>^"^ = i T^p.e^"" <?"^ ^'^T = i (trT) I, 

(2.9) (^yf"'!^ = i e^/^M j-i'^^^ (sj-pT ^ _ (2>ji 

(2.10) (i>r"^ = ^ T^^^^e'^)^'^ - i T^p^e^"'" 5"'', ^^^T'^ = ^^'T, tr (^'T = 0. 

2.2. The choice of the gauge. In the recent papers [31 13], the propagation of 
waves was studied in the structured medium under consideration. The results 
obtained describe the elastic waves in the two complementary pictures which arise 
due to the two different choices of the variables allowed by the gauge freedom (|1.12p - 
(|1.13p . In [3], the coframe gauge is used when the connection is zero (completely 
gauged out). The components of the coframe are the only dynamical variables. 
Geometrically, they represent an arbitrary orthogonal 3x3 matrix. Since the 
connection is trivial in this gauge, the torsion 2-form reduces to T" = d-d". Using 
the components of the coframe and frame, t?" = hfdx^, Ca = h^^di, we have for the 
holonomic torsion 

(2.11) T' = h^T"' = hl^d{h°') A dx^ = h^dkh'^ dx'' A dx^ . 
The dual 1-form then reads 

(2.12) r = Vj dx' = Kdkhf e'^^jdx^. 

Here e'^'j = h'^h^ph'Je"^ ^ is the Levi-Civita tensor with respect to the holonomic 
coordinate basis. 

A complementary approach is developed in [4] , where the only dynamical variable 
is the connection, whereas the coframe is fixed to its trivial constant value 1?" = 
(5" da;'. The torsion (|1.2p reduces in this gauge to T" = F^" A 1?^, and since in the 
Weitzenbock space the connection is given by (|1.14p . we have 

(2.13) r = SIT" = 5j.A", d{A-^yp dx' . 
Denoting := 5\^h°-a, we thus have 

(2.14) r = 1*^.^(1*^ A dx^ = uidkuf dx'' A dx\ 

This is completely equivalent to p. lip , with the difference that the orthogonal 
matrix /i^, representing the coframe, is replaced with the orthogonal matrix u := 
M^, representing the connection. 
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In [4], there are two technical deviations as compared to the model developed 
in [3]. The first deviation is the different choice of the variables mentioned above. 
Noticing that (u^duYkj = Uadku" is skew-symmetric in i,j (being an element of 
the Lie algebra of the orthogonal group), in [4] one chooses as a basic variable a 
new 3x3 matrix 

(2.15) A==*(u'^au), 

which has no a priori symmetries. This 7k- denotes the dualization operator which re- 
lates skew symmetric matrices to vectors, similar to the Hodge * operator. Writing 
out the indices explicitly, this matrix can be written in the following way 

(2.16) Ai.^^eu^u^d.u'^. 

This quantity is known in the polar elasticity theory under various different names: 
the Nye tensor [3D] , the wryness tensor [T3] , the second Cosserat deformation tensor 
[19j . the third order right micropolar curvature tensor |51j . or torsion-curvature 
tensor (note that this is somewhat a misnomer since it neither directly relates to 
torsion nor curvature in the sense of differential geometry). The inverse of (|2.16p 
yields the components of the torsion tensor p.l4p 

(2.17) uldku'} = Aike''j. 

Substituting this into (|2.14p . we can find a relation between the the second rank 
tensor T (which we can view as a matrix in R^) to the matrix A, so that 

(2.18) T = A-tr(A)I. 

The second deviation of [3] from [3] is a different structure of the Lagrangian. 
Namely, the Lagrangian of a rotationally elastic medium is assumed in [J] to be a 
function of orthogonal matrix only. 

2.3. Kinetic energy and total Lagrangian. The Lagrangian 4- form V = L dtAe 
is constructed by taking L ~ L{u, d^u, diu) as a quadratic function of the irreducible 
parts of j4;fe. In addition, the kinetic term is chosen to be 

(2.19) L^'^ = \ir{u^u) = i(a*<)(a,<). 

This form of the kinetic energy is well motivated, since when linearised it yields 
angular kinetic energy or rotation energy. Strictly speaking, in general the kinetic 
energy should be constructed by taking the microinertia into account |19[ 144] . 

Introducing, analogously to (|2.16p , the velocity of the deformations of the mate- 
rial continuum 

(2.20) Ait = \eu^v}^dtu^, 
and denoting ^4'* = g^^'Akt, we can recast (|2.19p as 

(2.21) L'^'" = i(ai<)(atO = A,,A'*. 

Taking into account the identity (|2.3p . the potential energy contains just two 
quadratic invariants. As a result, the general Lagrangian reads 

(2.22) L - L'^'^iAu) - LP°\Aik), LP°\Aik) = Xi (A^.)' + A2 A[ik]A^"'l 

The coupling constants Ai = |(c3 -f ^ci) and A2 ~ C1+C2 are constructed from the 
elastic moduli (|2.ip . 
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One can use different parametrizations of tlie ortliogonal matrices. For example, 
in [3], tfie spinor parametrization was used. Here we find it more convenient to 
describe an arbitrary orthogonal matrix with the help of the three real functions 
/?' which is based on the identification of S0{3) with the three-sphere. From these 
independent constituents, an orthogonal matrix is constructed as follows 



(2.23) 



{S^ + 2/3,(3' - 25^13^ + 2a(3''e'jk) 51, 



Here /?' = 5^-' Pj, + (3"^ ~ 1 (with — Sij(3'^(3^). The inverse matrix reads 

(2.24) u] = [si + 2(3^/3, - 261^" - 2a(3''ei,k) 5]. 
Substituting (P?^ and ^(TI^i into ^^U^i and (E^Ol), we find 

(2.25) Aik = 2{eH,/3'dk/3^ +(3idka-adk/3i), 

(2.26) Alt = 2{eHjl3'dtp' + pidta-adtPi) . 

For small /3, the above formulas can be linearised, so that Aik ~ ~2dkl3i, An « 
—2dtl3i, and the linearised model is described by the Lagrangian 

(2.27) L « 4/32 _ 4Ai(div;3)2 - 2A2(curl;3)2. 
Here as usual div/3 = diP', and (curl/3)i = 82(33 — ^3/32, etc. 

2.4. The nonlinear equations of motion. The complete nonlinear equations 
are more nontrivial. Denote the derivatives 



dLP° 



dAik 



dL 



kin 



dAit 



(2.28) H"" 
Then the field equations read 

(2.29) idtH^' - d,,H"^)P,j + 2{H''Qt,, ~ W'Qmj) = 0, 



where we introduced 

(2.30) P^J = 

(2.31) Qt^J 

(2.32) Qkij 



Sill 



{Sij(3i ~ Sii(3j) 
{5,j(3i - Sti(3j) 



dtP\ 



The equations ([2^ are valid for any Lagrangian L = L^''^{Aif) - LP°\Aik) 
with an arbitrary dependence on the variables Ait,Aik. However, for the specific 
choice of the quadratic Lagrangian (|2.2ip and (j2.22p . we have explicitly 

(2.33) W = 2A'\ = 2Ai A\5'^ + 2A2 A^'^^ 
The tensor (j2.30p is invcrtiblc and the inverse reads 

(2.34) [p-^y^ ^aS^^ - e^^''(3n. 



One can easily check that Pij{P ) 



Multiplying (|2.29p with this inverse. 



we obtain another convenient form of the fields equations: 
(2.35) dtW' - dkH'^ + 2{W*Gt/ - H^'^Gk/) = 0, 
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where we have introduced 

(2.36) Gtj' = Qtji{P-^ f ^e,''dt{api)+p'dtPj~p,dtl3\ 

(2.37) Gkj' = Qkji{P-'f ^e/'dkiaPi)+p'dkPj-PjdkP\ 

Note that both objects are antisymmetric in 

The resulting system of nonhnear differential equations (obtained after substi- 
tuting (11321), (P3U)) - (P3^ and (12211), (p:^ into ([123)) is quite nontrivial. It is 
possible, however, to find simple solutions under the additional assumptions. 

3. Spherically-symmetric solutions - the soliton 

Let us now look for the configurations with a 3-dimensional spherical symmetry. 
The corresponding ansatz reads 



(3.1) 



13' 



cos w, 



a ~ smw, 



where the scalar function w — w(t, r) depends on time and on the radial variable 

r = V Xix"^. 



Then (f2:25l) and (|2:26)) yield 



(3.2) Aik 

(3.3) Alt 



XlXk 



Hik ^ cos W + —j- 



w — sm w cos w 



5ik _ xiXk 



2 — w. 



Hereafter the dot and the prime denote the derivatives with respect to time and 
radius, respectively. 

Accordingly, the trace and the skew-symmetric parts of p.2p are 



(3.4) 
(3.5) 

As a result, we find 



sm w cos w 



A 



[Ik] 



2eikn cos w. 



(3.6) 
(3.7) 



4x' 



= 2Ai S''' 2w' - 4 



sm w cos w 



-4\^s'kn Xn _2 



cos w. 



It is straightforward to see that 
(3.8) 

Gkj 



4t' 
r 



tj 



0, 



jSmwcosui xiXk / sm If cos If 

fcj z ^ ~::2 



(3.9) 



-f- [x'Skj - XjSl) 



cos^ w 



Using these results, we find that the field equations (|2.35p reduce to 



(3.10) 



Ax' 
r 



?i — Ai I w" H — w' 
r 



■U{w) 
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where 

(3.11) U{w) = sm(2u;)[(A2 - Ai) + (As - 2Ai) cos(2w)]. 

Since the two terms in the round brackets arc the Laplacian in spherical coordinates, 
this equation of motion can also be written in the neat form id — XiAw = U{w)/r'^. 
Let us introduce a new function Lp := wr, and next let us rescale r i— > \/Xit, 
ip ^ \/AiV5, then p.lOp becomes 

U{^/r) 



(3.12) 



\3/2 

r 



= 0. 



The resulting equation is closely related to the spherical 'sine-Gordon' equation, 
see for instance [lOj [53] and references therein. It is worthwhile to mention that the 
model of Kontorova-Frenkel [37l [20l [5] was the first theory discussing media with 
dislocations where the 'sine-Gordon' equation and its soliton solutions naturally 
appeared. The main difference between previously studied spherical sine-Gordon 
equations and our equation p.l2p is that the nonlinearity carries an extra factor of 
1/r. This is similar to angular momentum when studying the spherical Schrodinger 
equation or Newton's equation. This additional factor has some interesting impli- 
cations. At large distances from the centre this term becomes negligible and asymp- 
totically we recover the wave equations. This analysis renders a strong support for 
the existence of soliton solutions in our model with a localised configuration near 
the centre. Below wc confirm this by numerical integration. 



w(r) 
l.O ^ 



0.4 I 



Figure 1 . Soliton solution of the rotational elasticity. 



In the static case, the equations of motion reduce to the single ODE 
(3.13) Xi{r^w')' + U{w) ^Q. 

For A2 = 2Ai this further reduces to [r^w')' + sin(2w) = 0, whereas when A2 



Ai 



one is left with (r^w')' — i sin(4i(;) = 0. The qualitative analysis of the equation 
p.l3p reveals the existence of static solutions that vanish at r = and approach 
asymptotically 7r/4 at infinity for r 00. 
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Numerical integration is straightforward. The form of the solution depends on 
the coupling constants and on the initial value of w'(0). However, the qualitative 
behaviour remains the same. As a specific example, Fig. [T] presents the soliton for 
A2 = Ai when w'{0) = 1. 

In the static case we can introduce the new function /, defined by 

(3.14) ^ ^ arctan{sinh(/(log(r)))} 

which transforms (j3.10p into an autonomous second order differential equation 

(3.15) /^^ + /a(l-tanh(/)/^)-2sinh(/)-4tanh(/)+2^(sinh(/)+tanh(/)) = 0, 

where /3 = log(r). This equation can now be analysed using standard techniques 
from ordinary differential equations or dynamical systems, and we find that this 
static system has three equilibrium points 



(3.16) f = 0, sinh/ = ± ^f"\^ , 

with eigenvalues (0, —1) in all three cases which in turn yields interesting (in)stability 
properties. 



4. Discussion 

The rotational elasticity model has many features similar to the model suggested 
by Skyrme [SH]. Although the Lagrangians are different, the dynamics looks quali- 
tatively the same. In particular, the remarkable feature of the rotational elasticity 
is the existence of solitons which are close relatives to the Skyrmions. The topo- 
logical nature of such configurations (which we demonstrate below) guarantees the 
stability of solutions obtained in the sense that the variations of the coframe and 
connection cannot change the value of the conserved charge. Moreover, the moduli 
space of these solutions is obviously compact. 

Recently it was noticed |54| that in the framework of the gauge approach in grav- 
ity and elasticity (that underlies the model under consideration, see Sec. II. 3p one 
can define an identically conserved current 3-form that gives rise to the topological 
charge that naturally classifies the field configurations. Such a 3-current "lives" on 
the 4-dimensional spacetime that is constructed as a foliation where the t = const 
slices coincide with the manifold M. Specialising to the case of the Weitzenbock 
geometry with the fiat curvature (|l.ip . this topological current reduces to 

(4.1) J'°P ^ KJ a DKp" + ^ KJ a Kp-' A K^". 

Here Ka^ = — is the contortion defined as the difference between the 
Riemann-Cartan and the Riemannian connections. This 3-form is identically con- 
served on the 4-spacetime, (iJ*°P = in view of (jl.ip . As a result, we can construct 
the topological charge 

(4.2) q=_1_JkJa a K,^. 

The integral is taken over the 3- manifold J\4 . The topological nature of the charge 
(j4.2|) and current (|4.ip is manifested in the fact that the corresponding conservation 
law does not depend on any dynamics of the fundamental fields and Ta^ [54] . 
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The integral (|4.2|) is taken over the whole 3-space and the result is a constant 
for configurations that we studied in the previous section when w vanishes at the 
origin and approaches constant value at infinity. Specializing to the connection 
gauge p.l4ll . we have Ta^ = (the coframe is "gauged away" : hf = 5f) and hence 



Fixing the asymptotic boundary conditions so that p — Sfu'j^ approaches unity 
near spatial infinity allows for the standard compactification [54] of the space M « 
S'^. We thus verify that the charge (|4.2p . (|4.3p is a degree ("winding number") of 
the maps A : S"^ — > 50(3) which are classified by the integer numbers that are the 
elements of the third homotopy group Z = 713(6*0(3)). By direct computation we 
can verify that Q = 1 for the soliton solution described above. 

Construction of a multi-soliton generalisation is an interesting problem. One can 
approach this along the lines proposed in |54] . As a first step, we can replace the 
ansatz p.ip by taking instead of (|2.23p the orthogonal matrix that is a product of 
N factors 



The dot denotes the usual matrix product. Such an ansatz would obviously generate 
a topological solution with the higher charge Q = N. The construction of the multi- 
soliton configurations is a more nontrivial problem which one could try to solve by 
using the instanton technique. 

The realization of the current theoretical model in the condensed matter systems 
is an interesting physical problem. The geometric constraint (|l.ip of the vanish- 
ing curvature rules out the disclinations. It was shown in |54j that dislocations 
are also not relevant to the point-like soliton solutions. The solitons of this type 
are mentioned by Unzicker [64], who calls them a "Shankar monopole" [58] that 
describes a point-like defect in a A-phase superfiuid Helium-3, [BS] HHl ISS] • As we 
mentioned already, our microrotational elasticity model belongs to the class of the 
so-called micropolar elasticity theories. The relevant physical continua describe, 
in particular, liquid crystals with rigid molecules, superfiuids, rigid suspensions, 
blood fluid with rigid cells, magnetic fluids, dust fluids etc. Eringen [19] (see also 
references therein) extensively studied micromorphic, microstretch and micropolar 
elasticity models with other defects beyond the dislocations and disclinations. We 
refer again to the paper of Randono and Hughes [54] where the point-like topolog- 
ical solution (called there a torsional monopole) is depicted and a possible physical 
manifestations of such structures are qualitatively related to certain nontrivial elec- 
tronic behaviour of condensed matter systems with strong spin-orbit coupling such 
as (l-|-3)-dimensional topological insulators. A careful quantitative analysis of such 
possibilities will be left for a future study. 
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Q 
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